UMN-TH-18 11/99 
TPI-MINN-99/36 
MADPH-99-1125 
hep-ph/9908311 
August 1999 



On 

Os ; Variations of the Neutralino Elastic Cross-Section with CP 

Os 

Violating Phases 

< 

[ Toby Falk 

Department of Physics, University of Wisconsin, Madison, WI 53706, USA 



m 

OO 
O 
OS 
Os 



O 
43 



X 



Andrew Ferstl and Keith A. Olive 

Theoretical Physics Institute, School of Physics and Astronomy, 
University of Minnesota, Minneapolis, MN 55455, USA 



<D , Abstract 



We analyze the neutralino-nucleus elastic cross-section in the MSSM, including contributions 
from CP-violating phases, using the four-fermi neutralino-quark interaction. Over a wide 
range of the MSSM parameter space we show the variations in the cross-sections due to these 
phases. We further concentrate on the regions which are consistent with constraints from 
the electric dipole moment of the electron, neutron and mercury atom. In the regions we 
examine in detail, we find suppressions by up to a factor of two, while enhancements in the 
cross-sections are no greater than ~ 10%. 



1 Introduction 



Supersymmetry provides a framework for constructing new models for physics beyond the 
standard model. The motivations for supersymmetry are well known, and include a stabi- 
lization of the mass hierarchy problem and unification of gauge couplings at the GUT scale. 
In the minimal supersymmetric extension of the standard model (MSSM) [TJ, conservation 
of R-parity, which can be related to baryon number, lepton number and spin, guarantees 
that the lightest supersymmetric partner (LSP) is stable and provides a candidate for the 
non-baryonic dark matter in the universe. Neutralinos, mixtures of neutral gauginos and 
Higgsinos, are the best motivated LSP candidate in the MSSM ||. 

For a given set of MSSM parameters, the exact identity of the LSP is determined, its 
annihilation cross-section and hence its relic density can be calculated, as can its interactions 
with ordinary matter. The latter is key in efforts to detect MSSM dark matter in laboratory 
experiments || ||. Indeed, over the years a great effort has been made to study the elastic 
cross-section of the LSP with nuclei in order to predict the possible detection rates in cryo- 
genic detectors. Recently, it has been shown 0, [| [7J that the CP violating phases in the 
MSSM may have significant effects on the detection rate of the neutralino. Here we continue 
our investigation of the neutralino-nucleon elastic cross-section over a wide range of MSSM 
parameters and show that these phases can give significant, variations to the cross-section. 
We also take into account the experimental electric dipole constraints of the electron, neu- 
tron and mercury atom, and show that variations in the cross-sections of order 0.4 - 1.1 
persist. 

CP violating phases can arise from several sources in the MSSM: the Higgs mixing mass 
((/,) in the superpotential, and the supersymmetry breaking parameters which includes the 
gaugino masses (Mi), the soft scalar masses, the bilinear term (Bli), and the trilinear terms 
(Af). Not all of these phases are physical ||. For our purposes, the phases associated 
with the gaugino masses are rotated away and Bli will be taken as real so that the vacuum 
expectation values of the Higgs fields are real. Also, the phases of the trilinear terms, Af 
(with the exception of At), will be set equal to one another, for simplicity. The trilinear stop 
mixing parameter A t has a (real) quasi-fixed point in its evolution from the GUT scale to the 
electroweak scale, so we take A t real. Then, the only physical phases in the MSSM are the 
phase associated with \x and the phase associated with the Af (6a)- Typically, in order 
to satisfy constraints]^] on the electric dipole moment (EDM) of the electron, neutron and 
199 Hg atom, CP violating phases must take on very small values (on the order of 1CT 2 — 1CT 3 
for MSSM parameters on the order of 100 GeV). This limit can be weakened ED, O if the 
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SUSY masses are taken to be upwards of a TeV, which is not a cosmological viability if the 
neutralino is gaugino-likell^]. Models with heavy first two generations of sfermions and a 



light third generation |T3|, on the other hand, may allow large phases while permitting the 
neutralino to annihilate into e.g. tau pairs, so that the cosmological constraints are satisfied. 
Alternatively, with two CP violating phases and several contributions to the electric dipole 
moments cancellations can occur and allow the CP violating phases to take on larger values 
TM - [113] . We neglect the (typically small) phase misalignment between the Higgs vevs, which 



is induced by one loop corrections to the Higgs potential in the presence of CP violation in 
the stop sector [19|, |2(J. The effect of this misalignment may be significant for large \/j 1 A t \/m\, 
particularly at large tan (3. 

Previous studies |5], ||] have shown that the CP-violating phases can have a large effect 
on neutralino-nucleus elastic scattering cross-sections. In certain regions of the parameter 
space, the elastic cross-section with matter can show a dramatic dependence (usually a 
minimum) as a function of the CP violating phase 9^. Here, this effect is studied in more 
detail. In particular, we determine the dependence of the cross-section on the phases across 
the M 2 -/x plane in the MSSM, including a separate investigation of the phase dependence of 
the spin-dependent and spin-independent cross-sections. We also examine the effect of the 
constraints of the electric dipole moments on the variation of the cross-section with respect to 
the phases. Generically, the EDM constraints require the phases to be sufficiently small that 
the cross-sections are only very mildly affected. Nevertheless, we do find parameter regions, 
predominantly where the neutralino is gaugino(bino)-like, where the constraints on the CP 
violating phases can be satisfied, and where the phases still produce significant changes to 
the scattering cross-section. In these regions, we find that the elastic cross-section may be 
enhanced by ~ 10% or reduced by as much as a factor of 2.5. 

In Section 2, the effective neutralino-quark interaction Lagrangian is presented and the 
dominant terms in the direct detection rate are discussed. In Section 3, the analytic ex- 
pression for the neutralino-nucleus scattering cross-section is found using the low energy 
effective Lagrangian. Finally, in Section 4 the variation of the elastic cross-section with the 
CP violating phases is shown. First, the effect of the phases on the spin-dependent, spin- 
independent, and total scattering rates is displayed for a particular set of parameters on the 
Mi — fi plane. For these plots, no EDM constraints are imposed so that the general behavior 
with respect to the CP violating phases can be shown. Then a general scan of parameter 
space is performed while imposing the electron, neutron and mercury EDM constraints so 
that a consistent and realistic understanding of the effect of CP-violation can be shown. 
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2 Effective Lagrangian 



The neutralino LSP is the lowest mass eigenstate of the linear combination of the bino, wino, 
and Higgsinos. In our notation, 

X = Z xl B + Z x2 W + Z X A + Z x4 H 2 (1) 

The neutralino mass matrix in the (B, W, Hi, if 2 ) basis is: 



N 



[Mi -M z sin W cos (3 M z sin W sin (3 \ 

M 2 Mz cos 6w cos (3 —Mz cos Ow sin f3 

—Mz sin 6w cos f3 Mz cos cos (3 — /i 

V M z sin 6^ sin /3 — M z cos 6^ sin (3 —/jl / 



where Mi and M 2 are the U(l) and SU(2) gaugino masses, tan (3 is the ratio of the Higgs 
vacuum expectation values, and fi is the Higgsino mixing mass. We assume a common 
gaugino mass at the GUT scale, which gives Mi = |tan 2 #iyM 2 , and we take tan/3,//, and 
M 2 as free parameters. The mass matrix is diagonalized by a matrix Z, Z*NZ~ 1 , which 
determines the eigenstate composition in ([!]). 
The sfermion mass 2 matrix can be written as 



2 , 2 2/3(T 3/ - Q f sin 2 6W)M| 

m 7 m f e~^f M'^ + mi + Qf sm A 9 W M' Z cos 2/5 , 



/ Ml + m) + cos2/3(T 3/ -Q/ sin 2 9 W )M 2 Z -m f m s e^f 



Here, Mum are the soft supersymmetry breaking sfermion masses which are real since we 
have assumed that they are generation independent and generation diagonal. Also, 

m f e^f = R ff i + A* f = R f \fi\ e w » + \A f \ e~ idA (2) 

where rrif is the mass of the fermion and Rf = cot/3(tan/3) for weak isospin +1/2 (-1/2) 
fermions. For simplicity, we will choose all Af, except for A t , to be degenerate at the weak 
scale, in our analysis. When the top quark Yukawa coupling is large at the GUT scale, the 
initial value of A t at M GUT tends to be damped as A t is evolved, and mi/2 sets the scale of 
A t at the weak scale. A t is therefore set equal to its quasi-fixed point value ~ 2mx/ 2 |[21|| . 
Here mu2 is the common gaugino mass at the GUT scale and is equal to ^M^olqutI 9 2 ■, 
where we've taken l/a GUT = 24.85. Finally, the sfermion mass 2 matrix is diagonalized by a 
matrix rj, •y]M 2 i]~ x , and can be parameterized by an angle Of with 



cos Of sin Ofe n f \ 




( Vn 


V12 \ 


sin Of e~ %l! cos Of J 




V V21 


V22 J 
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From the MSSM Lagrangian, we deduce the low energy effective four-fermi Lagrangian: 

L = xYl 5 XQa^u + «2i7 5 )<?i + UziXXMi + a4iXl 5 XQil 5 C[i + a 5 iXXQa 5 Qi + a&iXl^XMi (3) 

The Lagrangian is to be summed over the quark generations, and the subscript i refers to 
up-type quarks (i=l) and down-type quarks (i=2). The coefficients a are: 



OL\i = 



A{m 2 u - m 2 ) 



+ 



9 Z 



4m? cos 2 6L 



rv 1 2 I ry I 2 

^X3 I "~ 1^X4,1 



A(m 2 2i - m 2 x ) 
T 3i - 2ei sin 2 9 V 



(4) 



«2i = 



4(m 2 i - m 2 x ) 
9 2 

Am 2 cos 2 Q, h 



l^| 2 +l^| 2 ' 



+ 



J X4\ 



A{m 2 2i - m 2 ) 
2 



i^i 2 + iw,r 



(5) 



«3i = 



2(mfj — m 2 



-Re [(X,) {Ytf\ 



2(m 2 2i - m 2 ) 



grrigj 
Am m Bi 



Re (S u [gZ X 2 - g'Z xl \) Did 



+ 



1 



m 2 Hi m 2 H2 



+ Re (5 2i [gZ X 2 - g'Z xl ]) 



Df 



+ 



Cf 



™H 2 m %. 



(6) 



a 4i = 



2(mfj — m 2 ) 

gm qi Aj 
Am w B(m 2 Ht 



Re[(X t ) 



2(m|j - m 2 



-ite[(Wi) (ty* 



+ Ai?e (5 2i [gZ x2 - g'Z xl ])] 



(7) 



^ = -M /m,(yiKl!)V M /m[( " )W1 
+ a^^r ^ BlIm ^9Z X 2 - g'z xl \) 

- AJm(5 2l [gZ x2 -g'Z xl })} (8) 



«6i 



2(m 2 



li 



mi 



-Im[QQ (XiY 



2« 



mi 



-im[{vi) my 



igm q i 
Am w Bi 



Im (S u [gZ x2 - g'Z xl }) Did 



1 1 



m 



Hi 



m 



H 2 , 



+ Im (5 2i [gZ X 2 - g'Z xl ]) 



+ 



C? 



m 



Hz 



III 



Hi, 



(9) 



Where 



Yi 
Wi 
Vi 



* g m Qi z 



2m w Bi 



~ Vl2 e i9 Z 



* (Vi try rj, ry \ . * 9^ Z X 5~ 

gm q .Z* 5 -i 
V21- 



2m w Bi 



* f rz* 

V22 e i9 Z 



V22 



gm qi Z 



\5-i 



2m w Bi 



+ V21 [ i^g' z xi + g T 3i z X 2 



(10) 



Here, Sn is Z x3 (Z x4 ), 5 2 i is Z xA (— Z x3 ), £>j is sin/5 (cos/5), Aj is cos/3 (— sin/5), Cj is 
sin a (cos a), and .Dj is cos a (—sin a) for up (down) type quarks. The masses, mi and m 2 
correspond to the two squark mass eigenstates. m Hl 2 3 are the two scalar and pseudoscalar 
Higgs masses respectively (H 2 is the lighter scalar), a is the Higgs mixing angle. We note that 
(§) corrects an error in ||, where the expressions for a 2 and a 3 were presented previously. 
Hi is the hyper charge defined by = T 3i + ?/i/2. The expressions for the «j agree with 
those in ||.[] Also, in the limit of vanishing CP- violating phases, the spin-independent term, 
proportional to a 3i , and the spin-dependent term, proportional to a 2i , agree with those in 
I and H. 

From the Lagrangian we can calculate the cross-section of a neutralino scattering off of 
a target nucleus. All the terms in the Lagrangian, however, do not contribute equally to the 
cross-section. For example, 



\/f\ - - |-\|2 I > |2|| |2 (~Pxf ■ Pxi + m l) i-PNf ■ PNi + m N) 



mi 



ni 



AT 



\lt\ - - I -X 1 3 I /■ |2|, i 2 ( — PJV> 'PNi +™>lf) 



n? 



(12) 



JV 



x Our expressions agree with the corrected expressions in a revised version of || as communicated to us 
by the authors. 
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2 \/- |2| , 2 (-Px/ 'P* + m 



(/l«6iX7 X3tftl«) °c |Ci| law | i — (13) 



m 



where the pj are the four-momenta of the particle %. In the non-relativistic limit, the expres- 
sions (O)-(O) have a leading term proportional to the three momentum of the neutralino 
squared, the three momentum of the nucleus squared, or both. Here, Q = (N\q~iqi\N) , and 
($i = (N\qij5qi\N) , which corresponds to the quark contribution to the nucleon spin and 



is discussed in detail in pi. The axial- vector term represented by an is also found to be 



negligible 



K/|aiiX7^75X%7^K)| 2 oc 

, , 2 PxfPNf (PxiPm) + PxfPN t (PxiPN f ) + PN f PN z m 2 x - P Xf P Xi m% - 2m 2 x m 2 N 

Qfli n — 2 l i4 J 

m x m N 

In the non-relativistic limit, this reduces to a combination of the three momentum squared 
of the neutralino and nucleus. These terms are, hence, negligible compared to the terms 
multiplied by and which are not directly proportional to any three-momenta in the 
non-relativistic limit . 



The expressions for a 2 i and a^i in equations (|5|) and @ therefore, dominate the elastic 
cross-section. These expressions receive contributions arising from squark, Z°, and both 
scalar Higgs boson exchanges. The term proportional to «2j (equation (^j)) is a spin depen- 
dent coupling of the LSP to the nucleus due to the 7^75 factor, which is the spin projection 



in the non relativistic limit p4 |. On the other hand, the term proportional to (equation 
(§)) results in a coherent scattering of the LSP with the nucleus and as such is proportional 
to the mass of the nucleus. For large A nuclei, this term may be large for much of the MSSM 
parameter space. The steps to finding the cross-section from the effective Lagrangian are 
laid out in 0] and outlined below. 

3 Evaluating the Cross-Section 

The spin dependent cross-section can be written as 

32 



a 2 = —G 2 f m 2 A 2 J(J + 1) (15) 

7T 3 



where m r is the reduced neutralino-nucleus mass, J is the spin of the nucleus and 



A= j(a p (S p ) + a n (S n )) (16) 
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where 

a2i a GO „ _ a2i a W 



P i V2G f 1 i V2G f 1 V ; 

The factors A J - p,n ' ) depend on the spin content of the nucleon and are taken here to be 
Af } = 0.77,-0.38,-0.09 for u,d,s respectively and = A$\aP = A^,A^ = A® 
p6| The expectation values of the spin content of the nucleus, (S p>n ), depend on the target 
nucleus. Our results are for 73 Ge, which has {S p ^ n ) = 0.011,0.491, and for 19 F, which has 
{Sp,n) = 0.415, —0.047. The reader may refer to [|J] for details on these quantities. 
The spin-independent cross-section can be written as 

a 3 = ^[Zf p + (A-Z)f n } 2 (18) 

7T 

where 

— = E /t^ + I^E- (19) 
777- j q m a 27 ^ , m a 

P q=u,d,s 1 c,b,t 1 

and f n has a similar expression. The parameters fj?q are defined by {p\m q qq\p) = m p fj? q \ 
while ffy = 1 - E q = u ,d,sfTq WH- 0ur results used f!$ = 0.019,0.041,0.14 for u,d,s and 
/|J = 0.023,0.034,0.14 ||. There are other contributions to the spin-independent cross- 
section that depend on loop effects for heavy quarks and twist-2 operators which have been 
shown to be numerically small [[^| and are not included here. Note from (Q) and (fl0|) that 
the spin-independent cross-section is suppressed by the small nucleon mass. However, in 
contrast to the spin-dependent cross- sect ion, the spin-independent cross-section scales with 
the atomic weight of the scattering nucleus and therefore can dominate for heavy nuclei. We 
show results for 73 Ge, which has predominantly spin-independent interactions, and 19 F, for 
which either interaction can dominate, depending on the parameter region. 

The CP-violating phases can impact the cross-section in several ways. The phases can 
change the mass of the LSP. They will also change Z xi , which indicate the identity of the LSP, 
as well as the rjij, which indicate the mass eigenstates (and eigenvalues) of the sfermions. To 
get an understanding of how these effects arise, consider the case of the LSP being a nearly 
pure bino. (The neutralino can account for most of dark matter when it is mostly a bino.) 
In this case, Z x \ is nearly 1 and the other Z xi are very small. Then, the spin dependent term 
is approximately 



4(mfj — m y ) 



4(m| - m\\ 



* Vi I rv 

Vn^gZ x i 



* Vi I ry 



* f rv 

Vl2 e i9 Z 



+ 



Wi9 Z xl 



(20) 



7 



Hence, in this limit, the spin dependent cross-section will depend on the phases primarily 
through their effect on the mass of the neutralino, although the magnitude of the rjij will 
vary significantly, as well, if Ri\fi\ is comparable to \Ai\ (see (0)). 

To approximate the spin-independent term is somewhat more complicated due to the 
light scalar Higgs exchange. Keeping the leading terms we find 



«3i 



2« 



-Re 



( * try* \ ( * IJi I ry 

[-r] 12 e ig Z xl ) [riu-gZ. 



+ 



2 (ml; - ml 

gg'mqi 

Am w m 2 H2 Bi 



-Re 



( * / ry* \ I * Hi I ry 

[-r] 22 eig Z xl J [V2i^g Z x i 



Re (S u Z xl ) Did + Re (5 2i Z xl ) D\ 



(21) 



Due to the relative size of the Higgs mass to the squark masses, and the smallness of r\\2 
to r]ii, the Higgs exchange term can dominate in a^. In this limit, the spin- independent 
cross-section will depend on the phases through the mass of the neutralino, the phases in 
Z x i, and the phases in Z x3 ^ which exhibit a strong dependence on 6^. It would be expected, 
then, that the spin-independent term should be more sensitive to the changes in the CP 
phases as we will show below. We have not included the effect of the CP violating phases on 



the light Higgs mass |[20|| , which will typically be shifted by a few GeV for the masses and 
phases we study. 



4 Results 

Using the analytic expressions for the cross-section, we first show a set of contour plots on 
the M 2 — n plane for specific values of Af, M L ( R ), tan/3, and 6a- In Figures 1 - 6, we display 
the total cross-section, spin- dependent cross-section, and spin-independent cross-section for 
both 73 Ge (Figs. 1-3) and 19 F (Figs. 4-6). In these plots, the electric dipole moment 
constraints of the electron, neutron and mercury atom have not been imposed in order to 
show the general behavior of the cross-section as 9^ changes. The second set of plots (Figs. 
7 - 10) are scatter plots over the MSSM parameter space projected onto the M 2 — \i plane, 
concentrating on regions where the electric dipole constraints are satisfied. These plots serve 
as an existence proof that new CP violating phases can have a significant effect on the 
scattering rates while still satisfying the EDM constraints. Again, separate scatter plots will 
be made for 73 Ge and 19 F. 
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4.1 General Behavior 



In Figs. 1 - 6, we set A f = 3000 GeV, tan/3 = 3, m = 100 GeV, m A = 300 GeV and 
8a = tt/2. A t is set to its quasi-fixed point value ~ 2m 1 / 2 . The value of A j is chosen 
for easy comparison with Figures 7-8, where large Af is necessary to satisfy the electric 
dipole moments constraints. Figures la-c show how the total cross-section for scattering 
off of 73 Ge changes as 9^ changes from to 7r/4. The most significant change occurs in 
the lower right quadrant. To determine what causes this change, the spin-independent and 
spin-dependent cross-sections are plotted separately in Figs. 2 and 3. Direct detection 
experiments typically are sensitive to primarily one or the other type of interaction, and it 
is important to understand how the phases affect the different cross-sections. In Figs. 2a-c, 
we see that the spin-dependent cross- section varies only very mildly with 9^. In contrast, 
Figs. 3a-c show that the spin-independent cross-section (which has a maximum near the line 
fi = M2/2, where the bino-Higgsino mixing is large) for 73 Ge does change significantly, by 
up to a factor of 2, as # M is varied from to 7r/4, and the effect is largest in the lower right 
quadrant of the plane. The total cross-section's change is due to the large change in the 
spin-independent cross-section. 

Figures 4-6 are the corresponding pictures for 19 F. Again, the spin- dependent cross- 
section does not change much as 9^ increases, while the spin-independent term undergoes a 
significant change. Since 19 F is so light, the spin-independent term is comparable to the spin- 
dependent term over much of the parameter space. The total cross-section is affected by the 
change in 9^ in those regions where the spin-independent scattering cross-section dominates 
(for large \fj,\ and small M 2 ) and only little affected in the regions where the spin-dependent 
scattering cross-section dominates (for large M 2 and small | /x | ) . 



4.2 Behavior with EDM Constraints 

When the CP-violating phases are non-zero, the MSSM predicts a non-zero electric dipole 
moment for the electron, neutron and the mercury atom. In fact, over much of the parameter 
space, the values chosen in the previous plots will not satisfy the experimental bounds for 
the electric dipole moments (< 10~ 25 e cm |K| for the neutron) , (< 4.2 x 10 _27 e cm |3l| for 



the electron), and (< 9 x 10 _28 e cm |32| for the mercury atom). 

Within the MSSM, there are several contributions to the electric dipole moment. For 
the leptons and quarks, electric dipole moments are induced at one loop via exchange of 
sfermions, charginos, neutralinos, or gluinos (for quarks). For the neutron, there are ad- 
ditional contributions to the electric dipole moment, beyond the contribution by the quark 
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EDMs, which include a gluonic contribution and a quark color dipole contribution [11], |33], [15 



For a given set of values of the CP-violating phases, there may be significant cancellations 



between the individual contributions [|14| , jig , jig] . In general, the electric dipole moment of 
the electron provides a more stringent constraint over the M 2 — [i plane than does the neu- 
tron electric dipole moment, especially when uncertainties in the computation of d n are taken 
into account ||34j| . It is however, important to have more than one independent experimental 
constraint on edms, in order to obtain limits on the two CP violating phases. In [J34|| , it was 
argued that the edm of 199 Hg (induced through T- violating nuclear forces [E5[]) could be used 
to constrain the MSSM phases. Although it may require some degree of fine tuning, for any 
given set of tan/5, fi, M 2 , m , 6U i; and 9^ it is always possible to satisfy both of the electric 
dipole moments by choosing large enough magnitudes for the trilinear mass parameters, Aj 
|36| , |16| ) p7| (more correctly, acceptable solutions exist for either 9 a or Oa —* 9a + n)- In the 
absence of tuning the sfermion masses, the phases must take on values of order 10~ 2 — 
in order to satisfy the EDM constraints. 

For figures 7 - 10, a general scan of parameter space was performed over the M 2 ,/^, A, 
and rriQ while setting A t again equal to the quasi-fixed point. We have again taken tan (3 = 3 
and rriA = 300 GeV. The range of parameters is chosen via several criteria. From LEP2's 



search for charginos, there is a lower bound on the mass of the chargino of 95 GeV p8 |. This 
implies a restriction on the lower values of M 2 and fi of around 100 GeV. We also insist that 
the neutralino be the LSP, and we employ the recent LEP2 constraint on sfermion masses 
and discard runs with mj < 80 GeV |38| . Furthermore, when the sfermions are very massive 
and the LSP is gaugino-like, the neutralino's relic density is too large to be compatible with 
the observational lower limits on the lifetime of the universe [09|1 . Therefore for simplicity, 
we restrict the values of m-o so as to give cosmologically reasonable relic densities for all 
neutralino compositions (although of course mo is not constrained for Higgsino-like LSP's). 

In each figure, 9a and 9^ are fixed, and a scan is performed over M 2 , A, and m . If a 
set of parameters satisfies the EDM constraints, the value of the cross-section is calculated 
and compared to the value of the cross-section when 9a and 9^ are both zero. This scan is 
then projected onto the M 2 — /i plane. Figs. 7a and 7b are for 73 Ge for # M = n/8, 9a = 37r/8 
and 9^ = tt/4,9a = vr/2, respectively, and 8a and 8b are the corresponding figures for 19 F. 
In Figures 7 and 8, we require both the electron and neutron edm limits to be satisfied. For 
Figs. 7a and 8a, the range of parameters is: 200 GeV < M 2 < 800 GeV, 100 GeV < // < 
1000 GeV, 100 GeV < m < 200 GeV, and 2200 GeV < A < 3500 GeV. In the scan, M 2 
is incremented every 20 GeV, /i every 20 GeV, mo every 25 GeV, and A every 25 GeV. For 
figures 7b and 8b, the range of parameters is: 150 GeV < M 2 < 400 GeV, 100 GeV < fi < 
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1000 GeV, 100 GeV < m < 200 GeV, and 2700 GeV < A < 3500 GeV. In the scan, M 2 
was incremented every 10 GeV, fi every 20 GeV, mo every 10 GeV, and A every 25 GeV. 
We bin the resulting <t(# m , 9a ^ 0)/a(9 fJi , 9 a = 0) in steps of 0.1. A single symbol at a point 
means that every parameter set satisfying the EDM constraints has a cross-section ratio in 
the specified range. Some grid points have two symbols. Notice for all of the diagrams, there 
are regions where the cross-section is reduced. For some of the diagrams, however, there 
are regions where the cross-section is enhanced. In figures 7b and 8b, when 9^ = 7r/4 and 
Q a = tt/2 there are regions where the cross-section decreases by up to a factor of 2. Although 
this region of parameter space is small (in the region of 200 GeV< M 2 < 300 GeV and \i 
around 750 GeV), it has cosmological significance because it is in a region of parameter 
space where the neutralino is most nearly a bino, so that the neutralino can account for 
most of dark matter. By comparing the scatterplots to the previous contour plots, the 
largest changes identified by the scatterplots correspond to the region of parameter space 
that has the smallest cross-section for both Ge and 19 F. 

In |J, we had found cases in which the cross-section was reduced by over an order of 
magnitude for certain values of the phases. Those results show the maximum possible effect 
and such a reduction occurs at a different value of 9^ for each choice of M 2 and fi. Here, we 
have a milder effect, but one that should be taken as relatively generic. We have chosen two 
specific values of 6^, and scan the M 2 -[i plane. Therefore, in most cases, we will miss the 
value of 9^ for which the cancellation in the cross-section is at a maximum. Furthermore, 
we have lost a considerable amount of parameter space due to the imposition of the EDM 
constraints. This space could be enlarged if we were to relax the universality on the A-terms. 
By adjusting A e and A q separately, a larger portion of the plane would survive and larger 
effects would be visible. 

In Figures 9 and 10, we show the results when both the electron and mercury edms are 
satisfied. Figures 9a and 10a use the same range of parameters as those in figures 7a and 
8a. Likewise, figures 9b and 10b use the same parameters as those in figures 7b and 8b. The 
variation of the cross-section in these plots is very similar to the variation of the cross-section 
when the neutron and electron edms are imposed (figures 7 and 8) but a different section of 
the M 2 — plane is selected. These results (in figures 7 - 10) are in qualitative agreement 
with those of 0. 

It is important to note that we have chosen the ranges for our scans of A to coincide with 
the regions where cancellations are known to make the electric dipole moments small. As 
such, each point in Figs. 7 - 10a represents about 60 points in the scan (deprojected), and 
thus we find solutions for about 12% of the total of ~400,000 points in our scan. Similarly, 
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each point in Figs. 7 - 10b represents only about 10 points in the scan. In this case, we 
find solutions for about 1% of a total of ~400,000 points in the scan. These points serve 
to demonstrate the size of the uncertainties in the elastic scattering rates due to large CP 
violating phases. Again, the cosmological constraints may be satisfied by taking the first two 
sfermion generations heavy, while having a light third generation (13[. Lastly, we also note 
that there are regions at large mo(> ITeV) where the phases can be large without any tuning 
of the mass parameters. However, such large sfermion masses are cosmologically forbidden 
unless the LSP is mixed or Higgsino-like, in which case there tends to be an insufficient relic 
LSP abundance to make direct detection relevant. 



5 Summary 

We have made a comprehensive study of the effect of the two CP violating phases in the 
MSSM on the neutralino elastic scattering cross-section. Starting with the full 4-Fermi 
interaction lagrangian between neutralinos and quarks, we examined in detail the effect of 
the phases of the two dominant terms corresponding to the spin-dependent (0^2) and spin- 
independent (a 3 ) parts of the cross-section in the non-relativistic limit. 

We have shown explicitly the variation of the cross-section in the M2-/X plane for neu- 
tralino scattering off two sample nuclei, 73 Ge and 19 F. We showed the extent of the variation 
for the total cross-section as well as the individual spin-dependent and -independent contri- 
butions. It is the spin-independent piece in the bino-portion of the parameter plane that is 
most sensitive to the CP-violating phases. 

We have also considered in detail the regions of the parameter plane which are consis- 
tent with the experimental limits on the electric dipole moments of the electron, neutron, 
and mercury atom. The edms provide a very strong constraint and therefore, some degree 
of tuning is necessary to satisfy these constraints along with other cosmological and phe- 
nomenological constraints as well. In the surviving regions, the cross-section is typically 
reduced, by as much as a factor of ~ 2.5, while some regions in the M 2 -/i plane show slight 
enhancements in the cross-section by about 10%. 

It is also worth mentioning that the relic density of neutralinos is determined by their an- 
nihilation rate and the annihilation rate will also depend on the CP violating phases. In the 



general MSSM, these effect could be large [12] as the typical p-wave suppression of the anni- 



hilation cross-section is removed when the left and right handed sfermion masses are nearly 
degenerate. In the models such as the constrained MSSM (when universal soft masses at the 
GUT scale are assumed) this effect is small jnj. There is however, another contribution to 
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the annihilation cross-section which is unique to the presence of non- vanishing CP- violating 
phases. In the effective Lagrangian ([!]), terms which go as xi 5 X wn l De proportional to 
the three momentum of the neutralinos and will thus be negligible in the non-relativistic 
limit. However, the term proportional to will give non- negligible contributions to the 
annihilation rate for non- zero CP violating phases as can be seen by examination of eq.(|8]). 
In the non-relativistic limit we have 

E2 .2 of f \ 2 / \ 

-\a 5i \ m x \X - -J (22) 

where the total contribution should be summed over all fermions, f, which are less massive 
then the neutralino. Note that in the limit of vanishing CP violating phases, is zero since 
it depends on the imaginary contribution to (|I0D. A complete investigation of the effects of 
the phases on the total annihilation cross-section is left for future study. 
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Figure 1: The total cross-section for elastic scattering off of 73 Ge a.) for 9^ = 0, b) 9^ = ir/8, 
c.) 9^ = 7r/4. All plots have 9a = n/2. The axis are in units of GeV. The contours are in 
units o/HT 8 GeV~ 2 . 
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Figure 2: The spin- dependent cross-section for elastic scattering off of 73 Ge, for the same 
parameters as in Fig. 1. The contours are in units of 10~ 9 GeV' 2 . 
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Figure 3: The spin-independent cross-section for elastic scattering off of 73 Ge, for the same 
parameters as in Fig. 1. The contours are in units of 10~ 8 GeV~ 2 



18 




200 300 400 500 600 700 800 900 1000 200 300 400 500 600 700 800 900 1000 

|H| (GeV) W (GeV) 




Figure 4: The total cross-section for elastic scattering off of 19 F a.) for 9^ = 0, b) 9^ = 7i/8, 
c.) 6 ^ = 7r/4. All plots have 9a = n/2. The axis are in units of GeV. The contours are in 
units of KT 9 GeV~ 2 . 
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Figure 5: The spin- dependent cross-section for elastic scattering off of 19 F, for the same 
parameters as in Fig. 4- The contours are in units of 10~ 9 GeV~ 2 . 
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Figure 6: The spin-independent cross-section for elastic scattering off of 19 F, for the same 
parameters as in Fig. 4- The contours are in units of 1CT 10 GeV~ 2 . 
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Figure 7: The ratio of the elastic scattering cross-section with non-vanishing CP violating 
phases to the cross-section when the CP violating phases are zero, for 73 Ge. Here, the 
constraints from the eEDM and nEDM are taken into account and the scan of the parameter 
space described in the text is projected onto the M 2 — /i plane. In a) 9^ = it/8,0a = 37r/8 
and b) 6^ = tt/A, Q a = tt/2. 



22 



ct(9)/ct(9=0) 



a) 



800- 



700- 



600- 



§ 500- 



400- 



300- 



200- 



1 1 


"I"" 


t 




i'i'i'iii'i'i'i'iii'i'j; 


[1,1.1] 








===□=□===□=□====» 


□ □ 
[0.9,1] 








□ □□□□□□□□□□□□□□□□□□Q 


[0.8,0.9] 








□ □□□□□□□□□□□□□□□□□□□□□□□□a 






„□□ 


„„„ 


□ □□□□□□□□□□□□□□□□□□□□□□□□□X 




„„„„ 


„□□ 


„„„ 


□□□□□□□□□□□□□□□□□□□xxxxxxxx 










□ □□□□□□□□□□□□□□□□□xxxxxxxxx- 


+ + + + = o 


„„„„ 


„□□ 


„„„ 


□□□□□□□□□□□□□□□□xxxxxxxxxxx 
□ □□□□□ □□□□□□ □□□xxxxxxxxxxxx- 

□□□□□□□□□□□□xxxxxxxxxxxxxxx 










□ □□□□□□xxxxxxxxxxxxxxxxxxxx- 
□□□□□□xxxxxxxxxxxxxxxxxxxxx 










□□□□xxxxxxxxxxxxxxxxxxxxxxx 













100 200 300 400 500 600 700 800 900 1000 

|H| (GeV) 



ct(9)/ct(9=0) 



400- 



> 

& 300- 



iAeeeeeee 
^eeeeeee» 



b) 



i i i i i i i i i i 

[0.8,0.9] 

[0.7,0.8] 

e e 
[0.6,0.7] 

[0.5,0.6] 

AAAAAAAAAeeee&esee»T ttttii w ^ 

AAAAAAAAAgeg©eeeeeYTTT»T»TTY! 

xAAAAAAAAAeeeeeeeeTTTTTTTTTTTi 
AAAAAAAAAeeeeeeee»TT»T tttttt t> 

AAAAAAAAeeeeeeeeT»TYTTT»TTTYYl 



> y » e e s 

r e e e e e 



FQeeeeeAAA 



200- 



140 'i i i i i i i i i i' 

100 200 300 400 500 600 700 800 900 1000 

W (GeV) 



Figure 8: Same as Fig. 7 for 19 F. 
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Figure 9: Same as Fig. 7 using the constraints from the eEDM and the HgEDM. 
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Figure 10: Same as Fig. 9 for 19 F. 
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